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1 . 0  Introduction 

---  This  report  discusses  the  work  performed  in  the  second  year 
of  a  two  year  program  to  investigate  short-haul  non-li ne-of-sight 
optical  scatter  communication  links.  The  wavelength  region  con¬ 
sidered  for  this  study  is  the  middle  ultraviolet  (uv)  (2000  £  - 

o 

3000  A) .  This  region  is  attractive  for  communication  purposes 
because  the  absence  of  solar  background  radiation  makes  wide-angle 
receivers  practical,  furthermore,  the  short  wavelengths  result  in 
significant  scatter  from  molecular  atmospheric  constituents 
(Rayleigh  scatter) ,  so  that  the  availability  of  a  non-line-of- 
sight  link  is  extended  into  clear  weather  [1], 

The  goal  of  the  investigation  has  been  to  obtain  a  fundamental 
understanding  of  scatter  propagation  in  low  visibility  atmospheric 
environments.  In  particular,  the  work  has  focused  on  the  determina¬ 
tion  of  parameters  needed  for  optical  communication  system  design, 
such  as  angular  spread  and  multipath  time  spread. 

Since  low  visibility  environments  produce  significant  multiple 
scattering  effects  at  optical  wavelengths,  a  solely  theoretical 
approach  to  determining  these  propagation  parameters  is  extremely 
complex  [2,3].  Hence,  a  combined  experimental/theoretical  approach 
has  been  followed.  In  this  approach,  measurements  provide  order 
of  magnitude  ranges  for  the  propagation  parameters,  and  this  informa¬ 
tion  is  used  to  simplify  the  theoretical  propagation  equations. 

The  work  in  the  first  year  of  this  program  was  predominantly 
experimental  [4],  Experiments  in  Cambridge,  Massachusetts,  were 
conducted  to  obtain  representative  data  about  non-line-of-  sight 
propagation  of  middle  ultraviolet  light  in  city  environments. 
Experiments  were  also  conducted  to  measure  angu.l..r  spreading  beha¬ 
vior  on  a  line-of-sight  link.  Some  angular  spreading  data  was 
obtained  in  Cambridge,  but  the  majority  of  the  data  was  obtained 
during  a  field  trip  to  Lubec,  Maine.  In  Lubec,  the  frequency  of 
occurrence  of  maritime  fog  made  it  possible  to  collect  a  large 
amount  of  low  visibility  data  in  a  relatively  short  period  of  time. 

The  work  during  the  current  year  has  built  ^n  the  experimental 
data  base  of  the  first  year.  Experiments  were  performed  during  the 
current  year  to  supplement  this  data  base.  In  particular,  three 
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types  of  experiments  were  performed: 

1)  angular  spreading  experiments:  to  refine  the 
characterization  of  angular  spreading  begun  in  the  first 
year; 

2)  visible  vs.  uv  extinction  experiments:  to  measure  compar¬ 
ative  extinction  coefficients  in  fog  ?or  middle  uv 
light  (2537  A)  and  HeNe  laser  light  (6328  A)  . 

3)  multipath  time  spread  experiments  to  measure  pulse  distor¬ 
tion  through  fog  with  short  pulse  uv  sources; 

These  latter  experiments  will  not  be  reported  here,  since  they 
were  preliminary.  We  will  simply  note  that  for  optical  thicknesses 
less  than  10  on  a  0.43  km  path,  no  multipath  spread  was  observed 
with  a  2  usee  wide  pulse. 

In  addition  to  these  experiments,  work  is  near  completion  on  a 
modification  to  existing  measuring  equipment  which  will  allow  auto¬ 
mated  data  collection  for  future  angular  spreading  measurements  [5] . 

Theoretical  work  has  progressed  in  a  number  of  areas  during  the 
current  year.  Comparisons  have  been  made  betv/een  the  measured  data 
and  predictions  from  single  scatter,  diffusion  and  multiple  forward 
scatter  theories.  The  latter  theory,  which  is  on  appl ication  of 
results  from  turbulence  theory  to  a  scattering  environment  [6],  is  in 
reasonably  good  agreement  with  the  measured  data.  Calculations  from 
various  theories  have  been  supplemented  by  Monte  Carlo  simulations 
where  necessary. 

Some  progress  has  been  made  in  utiliz  ing  the  experimental  data 
to  simplify  the  linear  transport  equation,  which  describes  the  propa¬ 
gation  in  a  scattering/absorbing  medium  [7,8J.  As  initially  observed 
during  the  first  year  of  this  work  [4] ,  the  amount  of  light  which 
reaches  the  receivers  at  large  angles  from  the  line-of-sight  direction 
is  relatively  insensitive  to  random  spatial  variations  in  the  extinc¬ 
tion  coefficient.  This  fact  has  been  used  to  make  simplifications  to 
the  integral  equation  for  the  scattered  component  of  the  field. 

The  experimental  and  theoretical  work  conducted  during  this  year 
are  described  in  more  detail  in  the  sections  which  follow. 
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2 . 0  Experiments 

This  section  presents  the  results  of  the  angular  spreading  and 
visible  vs  uv  extinction  measurements.  The  presentation  of  the  data 
is  preceded  by  a  discussion  of  the  methodology  used  in  characterizing 
the  weather  conditions,  and  a  description  of  the  measuring  equipment 
utilized. 

2.1.0  Measurement  of  Weather  Conditions’ 

To  relate  theoretical  work  to  the  experimental  data,  it  is 
sufficient  to  know  the  scattering  coefficient,  absorption  coefficient 
and  the  single  scatter  phase  function  resulting  from  the  weather  con¬ 
ditions  prevailing  at  the  time  of  the  experiments.  It  is  a  relatively 
simple  matter  to  measure  the  scattering  and  absorption  coefficient, 
particularly  at  middle  uv  wavelengths  where  the  primary  absorber  is 
molecular  ozone.  If  ozone  concentration  is  measured  at  the  time  of 
the  experiment,  then  a  single  extinction  measurement  will  determine 
tne  scattering  coefficient. 

In  practical  situations  of  cours©  at  any  time  the  scattering 
coefficient  varies  randomly  as  a  function  of  distance  along  the  path. 
Thus,  an  extinction  measurement  involving  integration  over  the  entire 
path  is  necessary.  In  general,  a  path-  integraLed  ozone  absorption 
coefficient  would  also  be  necessary.  However,  if  the  ozone  concen¬ 
tration  in  a  region  is  relatively  stable  (e.g.,  there  are  no  ozone 
sources  nearby  )  ,  and  the  path  is  not  too  long,  a  simple  point 
measurement  of  ozone  concentration  suffices. 

Unfortunately,  there  is  no  easy  way  to  measure  the  single 
scatter  phase  function.  Measuring  the  phase  function  directly,  or 
measur ing  the  particle  size  distribution  -  from  which  the  phase 
function  can  be  determined  by  Mie  calculations  -  is  difficult  and 
time  consuming.  However,  if  observations  are  made  of  the  macroscopic 
features  of  the  weather  conditions  (such  as  qualitative  weather  type 
and  relative  humidity) ,  then  published  experimental  or  theoretical 
phase  functions  may  be  used  [9,  10]. 

The  procedure  used  to  characterize  the  atmosphere  during  our 
experiments  is  as  follows:  First,  a  qualitative  weather  type  (e.g., 
fog,  snow)  is  noted,  along  with  any  other  relevant  descriptive 
information  (e.g.,  high  variability  in  conditions).  Second, 
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temperature  and  relative  humidity  are  recorded  (sometimes  after  the 
fact) .  Third,  a  point  ozone  concentration  measurement  is  made  at 
the  receiver.  Finally,  the  average  path  extinction  coefficient  is 
determined. 

This  latter  quantity  is  derived  by  comparing  the  power  received 

on  a  clear  day  in  a  1  milliradian  (mr)  field  of  view  (FOV)  with  that 
received  in  the  same  FOV  under  low  visibility  conditions.  Since  a 
J.  mr  FOV  will  collect  predominantly  unscattered  light,  the  negative 
log  of  the  ratio  of  these  two  power  measurements  is  the  integrated 
optical  thickness,  t: 

p 

L  low  visibility 

x  =  /  a(z)dz  =  -In  — - -  (2.1-1) 

0  clear 

Here  a(z)  is  the  extinction  coefficient  as  a  function  of  z,  the 
distance  along  the  line-of-sight  path  from  source  to  receiver.  L  is 
the  pathlength.  t/L  is  the  average  extinction  coefficient  along  the 
path. 

2.2  Equipment  Description 

A  variety  of  source,  receiver  and  filter  combinations  were  used 
for  the  different  experiments.  Tables  2-1  through  2-3  list  their 
specifications.  Figure  2-1  shows  schematic  diagrams  of  the  three 
receivers.  Table  2-4  shows  which  combinations  of  equipment  were  used 
for  each  type  of  experiment. 

In  addition  to  the  equipment  listed  in  the  tables,  ozone  measure¬ 
ments  were  made  with  a  Mast  Development  Company  Model  724-5  Ozone 
Meter. 


Table  2-1  Source  Specifications 
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2.3  Description  of  Field  Trip  Experimental  Sites 

Almost  all  of  the  angular  spreading  data  was  obtained  during 
field  trips  to  the  Maine  coast.  In  the  first  year  of  the  program 
experiments  were  conducted  at  Lubec,  Maine  [4],  During  the  second 
year,  experiments  were  conducted  in  the  nearby  community  of 
Eastport,  Maine.  This  section  describes  the  experimental  sites 
for  the  Eastport  field  trip  (A  description  of  the  Lubec  site  is 
given  in  reference  [4].)  Section  2.3.1  discusses  the  prevailing 
weather  conditions  in  the  Eastport  area.  Section  2.3.2  discusses 
the  experimental  link  geometry. 

2.3.1.  Weather  conditions  in  the  Eastport  Area 

The  city  of  Eastport  is  located  in  Washington  County,  Maine, 
near  the  mouth  of  the  Bay  of  Fundy.  (See  Fig.  2-2).  This  area 
ranks  third  on  the  east  coast  in  frequency  of  fog  during  the  year, 
with  an  average  of  65  days  of  fog  [11] .  As  can  be  seen  from 
Fig  2-3,  the  highest  incidence  of  fog  is  during  the  summer  months, 
particularly  July  [12 ] . 

During  the  first  year  18  day  field  trip  to  Lubec  (July  24  - 
August  11,  1978),  we  experienced  7  days  of  heavy  fog.  Six  of  these 
days  were  consecutive  and  during  two  of  them  the  fog  did  not  lift 
for  the  entire  24  hours.  We  experienced  heavy  fog  during  9  days  of 
the  6  week  field  trip  to  Eastport  (June  9  -  July  20,  1979).  The 
relative  frequency  of  fog  in  Eastport  was  lower  than  in  Lubec  because 
the  Eastport  experimental  facilities  were  located  in  a  cove  partially 
sheltered  from  the  winds  which  blow  the  fog  inland. 

The  fog  type  encourtered  during  these  field  trips  is  classified 
as  a  stable  advection  fog  [13].  The  formation  mechanism  is 
typical  of  coastal  maritime  fogs,  in  which  warm,  moist  air  condenses 
over  colder  sea  water  and  is  then  blown  inland.  In  general,  the  fog 
occurinq  during  the  experiments  had  existed  for  some  time  before 
blowing  inland,  and  hence  was  a  "mature"  rather  than  "evolving" 
fog  ( 14 ] . 

The  ozone  concentration  measured  during  the  Eastport  field 
trip  was  between  0.5  and  2.0  parts  per  hundred  million  (pphm) ,  but 
the  concentration  was  usually  in  the  1-2  -  i.8  pphm  range.  The 
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Fiy.  2-2:  Location  of  Eastport,  Maine 
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Fig.  2-3:  Average  Number  of  Days  of 
Heavy  Fog 

(Visibility  <  0.5  miles) 
for  Period  Jan.  1960  - 
Dec.  1978 
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pathlengths  used  in  the  field  experiments  were  short  enough 
(<  1  km)  so  that  ozone  concentrations  in  the  above  ranges  did 
not  produce  significant  losses.  At  a  pressure  of  1  atmosphere 
the  ozone  absorption  coefficient  is  related  to  ozone  concentra¬ 
tion  by  [15] . 


aQ  (knf1)  =84.4 


(2.3-1) 


where  C~  is  the  ozone  concentration  in  pphm  and  T  is  the  ambient 
j  o 

temperature  in  K.  Table  2-5  shows  ozone  absorption  coefficients 


for  the  ranges  of  CQ  encountered  in  the  field  trips.  Also 


included  in  the  table  are  values  of  ozone  transmission  for  two 
representative  pathlengths  (0.43  and  0.88  km).  The  worst  case 
transmission  loss  is  less  than  a  factor  of  2.  This  is 


relatively  insignificant  by  comparison  with  the  scattering  losses 


of  102  -  105  in  fog. 


Table  2-5  Ozone  Absorption  Coefficients  and  Corresponding 
Transmission  Due  to  Ozone 

Cq3  (pphm) 
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2.3.2  Eastport  Link  Geometry 

A  detailed  view  of  the  two  Eastport  sites  is  shown  in  Fig  2-4 
and  the  topography  of  the  main  propagation  link  (at  the 
Washington  County  Vocational  Technical  Institute's  Marine  Center) is 
shown  in  Fig  2-5.  The  topography  of  the  Eastport  airport  (site  #2) 
is,  of  course,  flat,  with  a  rough  asphalt  surface.  For  the 
experiments  conducted  at  the  airport  site,  both  the  source  and 
receiver  were  approximately  5  feet  off  the  ground.  The  facilities 
at  the  Vocational  Institute  allowed  for  a  long  path  of  0.4  3  km*.  Two 
shorter  paths  were  also  established  by  mounting  sources  on  telephone 
poles  closer  to  the  receiver  location  at  the  end  of  the  dock.  The 
two  short  path  distances  are  .22  km  and  .32  km.  Experiments  at  the 
airport  were  conducted  with  a  fixed  source  and  the  portable  Honeywell 
receiver.  Pathlengths  varied  from  about  1.2  km  to  .88  km. 

2.4  Angular  Spreading  Experiments 

Although  a  significant  amount  of  fog  angular  spreading  data 
was  obtained  during  the  first  year  of  this  program,  most  of  it  was 
limited  to  a  relatively  narrow  range  of  optical  thickness  (5<t<9) . 
Thus  one  important  goal  of  this  year's  experiments  was  to  explore 
the  regions  x  <  5  and  t  >  9 . 

Another  major  goal  was  to  verify  a  number  of  basic  features  of 
the  fog  angular  spectrum  that  had  been  identified  after  analyzing 
the  Lubec  data.  Briefly,  these  features  are  [ 4 J : 

1)  There  is  a  measureable  unscattered  signal  component  for 
t  <  9; 

2)  This  unscattered  component  falls  below  the  scattered 
component  for  x  >  10  and  a  uniform  angular  spectrum 
result  (at  least  for  the  angular  region  j  y [  <  12  mr 

3)  For  t  >  5,  the  10  dB  hallwidth  of  the  angular  spectrum 
is  between  5  and  10  mr .  The  spectrum  tends  to  broaden  as 
x  increases; 

4)  For  off-axis  angles,  the  received  signal  is  relatively 
insensitive  to  spatial  variations  in  the  extinction 
coefficient,  while  the  on-axis  (or  unscattered)  signal 
is  quite  sensitive  to  these  variations 


Overhead  View  of  2  Eastp< 
Sites  (WCVTI  and  airport 


Fig.  2-5:  Side  View  of  WCVTI 
Propagation  Link 
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5)  The  relatively  uniform  portion  of  the  angular  spectrum 
extends  pver  a  considerable  range,  at  least  as  far  out 
as  25°. 

Pigs,  3-1  and  3-2  of  reference  [4]  are  pictorial  representations 
of  these  features. 

Finally,  it  was  our  goal  to  determine  to  what  extent  the 
above  features  were  influenced  by  absolute  pathlength.  Most  of 
the  measurements  to  date  had  been  at  a  single  pathlength;  optical 
thickness,  which  is  a  normalized  pathlength,  was  used  to 
characterize  the  path.  We  wanted  to  find  out  how  the  measured 
results  for  a  given  optical  thickness  would  be  affected  by  different 
pathlengrhs . 

Although  the  fog  was  never  thick  enough  during  the  Eastport 
field  trip  to  get  data  for  t  >  9,  most  of  the  other  goals  were 
achieved:  significant  data  was  obtained  for  x  <  5?  characteristics' 

1) ,  3)  and  4)  listed  above  were  repeatedly  observed  in  the  region 
5  <  x  <  9  (although  some  of  the  angular  spectrum  10  dB  half  widths 
were  less  than  5  mr) ;  some  multiple  pathlength  data  were  obtained. 

In  the  rest  of  this  section,  the  data  are  presented  and  discussed. 

Table  2.6  is  a  summary  of  the  small  angle  signal  versus 
collecting  FOV  data,  and  Table  2-7  summarizes  the  angular  scan  data. 
Fig  2-6  shows  a  comparison  of  signal  collected  by  a  narrow  FOV 
(1  mr  full  angle)  receiver  with  that  collected  by  a  wide  FOV  (84° 
full  angle)  receiver. 

In  order  to  see  how  the  angular  spectrum  broadens  as  t 
increases,  the  small  angle  signal  vs.  FOV  data  has  been  plotted 
with  x  as  a  parameter.  The  curves  are  drawn  throuyn  the  mean 
values  of  the  data  of  each  FOV.  (See  Fig.  2-7)  .  Observe  that  most 
of  the  data  are  consistent  with  the  conclusion  that  the  angular 
spectrum  broadens  as  the  optical  thickness  increases. 


ARO-18 


Table  2-6 

SUMMARY  OF  UV  SIGNAL  VS.  FOV  DATA:  EASTPORT,  1979  *^ 


Optical 

Multiplication 

of  Mean  Signal 

Date 

Time 

Thickness 

5  mr 

10  nr 

15  mr 

6/23 

4:10  “m 

3. 1-4. 4 

1. 7-2.1 

2. 4-3.0 

2. 8-3. 6 

4:00 

2.8 

2.1 

2.6 

2.9 

4 :  2C 

1. 8-2.1 

1. 8-2.4 

1. 8-2.4 

1.9-2. 5 

4:25 

3.4-3. 7 

1.8-2. 5 

2. 1-2.9 

2. 1-2. 9 

4:20 

5.9  (?) 

1.9 

1.8 

2.7-3. 5 

4:30 

4.0 

2.4 

3.5 

4.2 

4:33 

4 . 0 

2.0 

3.0 

3.4 

4:36-40 

4.3 

2.0 

3.7 

3.9 

5:15 

3.9 

2.3 

3.0 

3.5 

5:10 

3.4 

2.1 

2.8 

3.1 

5:20 

3.7 

2.5 

2.8 

3.1 

5:20 

4.0 

2.1 

3.2 

3.7 

10:20  pm 

6,3 

2.8 

4.9 

7.6 

11:40 

5. 6-6.1 

2. 4-3.9 

3. 4-5.6 

3. 6-5. 9 

11:45 

6.2 

3.1 

5.2 

5.9 

7/1 

12:20  am 

*  ★  * 

12:23 

5.2 

2.4 

4.0 

4.9 

12:30 

4.9 

2.1 

3.5 

4.0 

3:30 

4.8 

2.9 

5.6 

6.2 

3:40 

4. 4-4. 7 

2. 8-4.1 

5. 6-8. 2 

5. 5-8.0 

3:50 

4.6 

3.7 

5.8 

7.3 

11:10-15 

pm  3.3 

2.8 

3.9 

4.5 

7/2 

5:43  am 

6.5 

7.7 

12.5 

14.8 

5:45-48 

6.3 

3.8 

8.3 

11.7 

10:13  pm 

4.3 

1.8 

3.2 

3.9 

11:50 

*** 

11:53 

6.4 

2.3 

3.3 

4.9 

11:55 

6.8 

3.0 

4.8 

5.8 

7/3 

12:00  am 

7.1 

2.7 

4.0 

6.2 

1:35 

it  it  ★ 

7/16 

3:50  am 

4.0-4. 4 

1. 3-1.8 

2. 1-3.0 

— 

***Data  obtained  Dut  weather  too  variable  to  make 
,  a  reliable  estimate  of  optical  thickness. 

w  "f  Pathlength  *  0.43  km,  FOV  for  unscattered  measurement  =  1.0  mr. 


Table  2-7 


SUMMARY  OF  ANGULAR  SPECTRUM  DATA:  EASTPORT,  1979  * 


Date 

Time 

6/23 

10:25  pm 

10:50 

11:50 

6/24 

12:00  am 

6/25 

10:30  pm 
11:10 

7/1 

3:52  am 
4:15 

7/3 

12:07  am 
12:30 
12:45 
1:25 

7/12 

1:30  am 
1:50 

3:20 

3:50 

Optical 

Thickness 

Type 

5. 5-5. 9 

az 

8. 4-8. 7 

el 

.  5.3 

el 

5.1 

el 

0 

az 

0 

el 

4. 4-5.1 

el 

*** 

7.1 

az 

6.3 

el 

5. 1-6.0 

az 

5. 0-6.0 

az 

*** 

az 

6. 0-6. 5 

az 

7. 0-8.0 

el 

8.3 

el 

10  dB  Half  Width 

4  mr 

4- 6  mr 
3  mr 

2.5  mr 

1  mr 
1.25  mr 

5- 7  mr 


4. 5-5.0  mr 
7  mr 

(5  mr  FOV)*r - 

(5  mr  FOV)  - 

(5  mr  FOV) 

(5  mr  FOV)  — 

(5  mr  FOV) 

9  mr 


***  Data  obtained,  but  weather  too  variable  to  make 
a  reliable  estimate  of  optical  thickness. 

**  Scan  was  done  with  a  5  mr  FOV  instead  of  a  1  mr  FOV 


*  Pathlength  =  0.43  km 


Fiy.  2-6  Wide  FOV  (84  )  vs.  Narrow 
FOV  (1  milliradian)  Signal 
(Pathlength  =  0.43  Xm) 
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The  inversion  of  the  curves  for  6.0  <  x  <6.3  and  6.8  <  x  <  7.1 
appears  to  violate  this  conclusion,  but  the  large  error  bars  on  the 
6.0  <  x  <  6.3  data  overlap  (or  nearly  overlap)  the  error  bars  for 
the  6.8  <  x  <  7.1  data.  Indeed,  if  the  measurements  on  7/2/79 
(5:45-48  am)  are  not  included  in  the  data  for  the  latter  range,  the 
two  curves  are  quite  close  together.  .  There  are  other  examples  in 
Table  2-6  of  data  that  are  not  consistent  with  this  conclusion,  but 
the  majority  of  the  data  tend  to  support  it. 

The  curves  in  Fig.  2-7  only  show  signal  level  relative  to  the 
1  mr  value  in  any  optical  thickness  region.  Of  course,  'the  signal 
levels  for  the  larger  ranges  of  optical  thickness  are  below  those 
for  the  smaller  ones.  Figures  2-8  through  2-10  show  this  optical 
thickness  dependence  explicitly,  with  FOV  as  a  parameter  (The  5  mr, 
10  mr  and  15  mr  curves  are  shown  in  separate  figures  for  clarity.) 
Included  in  these  figures  is  the  corresponding  curve  for  an  84° 
full  angle  FOV,  obtained  from  Fig.  2-6. 

As  can  be  seen  from  the  figures,  the  data  points  for  any  single 
FOV  very  nearly  lie  on  a  straight  line.  The  slope  of  this  line  is 
dependent  upon  the  FOV,  ranging  from  3.8  dB  transmission  loss  per 
optical  thickness  (for  FOV  =  5  mr)  to  1.25  dB  per  optical  thickness 
(for  FOV  =  84°) .  Thus  the  data  suggest  a  general  exponential 
extinction  law  valid  for  a  receiver  processing  both  the  scattered 
and  the  unscattered  signal  components: 


P  -f(F°VU  (2.4-i, 

r  T 

where  f(FOV)  is  always  less  than  1  and  is  a  monotonically  decreasing 
function  of  the  FOV. 

We  defer  any  further  discussion  of  the  function  f (FOV)  until 
after  the  multiple  pathlength  data  is  presented.  Some  understanding 
of  the  pathlength  dependence  of  the  angular  spreading  is  needed  in 
order  to  use  data  taken  over  different  pathlengths  to  compute  a 
universal  function  f(FOV). 

Multiple  path  transmission  experiments  were  conducted  at  the 
Eastport  airport,  with  the  source  location  fixed  and  the  receiver 
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location  variable.  The  portable  Honeywell  sensor  was  used  with  both 
the.  22°  and  84°  field  of  view.  Measurements  were  made  with  both 
source  and  receivers  mounted  about  5  feet  from  the  ground. 

Time  limitations  due  to  weather  variability  made  it  impossible 
to  measure  optical  thickness  directly  during  these  experiments. 

(Aiming  to  within  a  milliradian  requires  a  few  minutes  of  time 
under  the  best  of  circumstances) .  Thus,  the  values  of  x  (and  hence 
a,  the  extinction  coefficient)  were  inferred  from  the  84°  FOV  trans¬ 
mission  data  by  comparing  the  data  at  a  pathlength  of  0.47  km  with 
the  narrow  versus  wide  FOV  data  in  Fig  2-6.  (The  comparison  here  is 
not  perfect,  since  the  latter  data  were  taken  over  a  0.43  km  path 
with  somewhat  different  boundary  conditions.  (See  Fig  2  -5b.) 

However,  the  0.04  km  difference  is  not  very  significant,  and  in  both 
sets  of  boundary  conditions,  most  of  the  scattered  signal  will  come 
from  the  upper  hemispheres.) 

Table  2-8  presents  the  results  of  the  multiple  path  transmission 
experiments  for  two  fog  thicknesses,  a  -  11.0  km  ^  and 

a  *  15.2  km  ^ .  In  Figures  2-11  and  2-12,  these  data  are  compared 
with  data  taken  at  the  same  fields  of  view  but  with  a  fixed  pathlength. 
For  the  84°  FOV,  the  fixed  pathlength  measurements  were  made  on  the 
0.43  km  path  in  Eastport.  For  the  22°  FOV,  the  fixed  pathlength 
measurements  were  made  on  the  0.30  km  path  over  water  at  Lubec. 

The  data  in  Figures  2-11  and  2-12  are  somewhat  equivocal  with 
respect  to  determining  the  sensitivity  of  the  result  to  pathlength. 

If  the  t  =  9.7  data  points  can  be  ignored,  then  the  remaining 
multiple  pathlength  data  lie  f.  rly  close  to  the  fixed  pathlength  line. 
However,  there  are  not  enough  data  at  the  larger  optical  thicknesses 
to  justify  calling  the  x  =  9.7  point  anomalous.  In  fact,  if  tha 
x  =  13.4  data  points  are  ignored,  then  a  very  different  conclusion 
results . 

There  is  further  evidence  of  the  relative  pathlength  insensitivity 
of  our  results,  however.  During  the  Lubec  field  trip  in  1978.  narrow 
versus  wide  FOV  measurements  were  made  with  a  17°  full  angle  FOV,  at 
both  0.30  km  and  1.6  km.  (See  Figures  A-18  and  A-19  of  reference 
[4].)  These  two  sets  cf  data  are  plotted  together  in  Fig.  2-13. 
Although  most  of  the  0.3  km  data  are  for  x  <  6  and  most  of  the 
1.6  km  data  are  for  x  >6,  it  is  clear  from  the  figure  that  both 
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TABLE  2- 8a.  MULTIPLE  PATH  TRANSMISSION  LOf,S  DATA 
(a  =  II. 0  km"1) 


Distance 

(km) 

T 

I84°  MB) 

X 22 ° (dB) 

*-ext(dB> 

.29 

3.2 

3 

5 

14 

.47 

5.1 

7 

10 

22 

.88 

9.7 

21 

23 

42 

TABLE  2- 8b. 

MULTIPLE  PATH  TRAN 

EMISSION  LOSS 

DATA 

(a  = 

15.2  km_1) 

Distance 

(km) 

t 

*Cq4°  (^B) 

0^22° 

*ext(dB) 

,29 

4.4 

5 

7 

19 

.47 

7.1 

10 

13 

31 

.88 

13.4 

20 

24 

58 

PRO- 2 3 
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V  SLOPE  s  1.5  dB/r 

Straight  line  fit 
to  data  in  Fig. A-20 
of  reference  [4]. 

\  (Solid  line  is 
v  region  actually 
N  measured. ) 

I 

I _ I _ L_ 

10  12  14 


Fig.  2-12:  Comparison  of  Multiple 
Pathlength  ,  Data  with 
Straight  Line  Fit  to  Fixed 

Pathlength  Data  (FOV  =  22° 
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groups  of  data  fall  along  a  straight  line,  whose  slope  is 
approximately  2.5  dB  per  optical  thickness. 

Fig.  2-13,  combined  with  Figs.  2-11  and  2-12,  is  strong 
evidence  that,  at  least  over  the  pathlength  ranges  covered  in 
our  experiments,  the  character  of  the  transmission  curve  is  not 
highly  sensitive  to  absolute  pathlength.  Table  2-9  summarizes  the 
available  data  for  transmission  vs.  FOV  in  terms  of  the  slope  of 
the  transmission  curve.  The  corresponding  value  of  f(FOV),  the 
coefficient  in  the  exponential  transmission  law  Eq.  2.4-1,  is  shown 
for  each  EOV. 

Figure  2-14  is  a  plot  of  f(FOV)  versus  FOV.  (Note  that  the  1  mr 
point  is  assumed  to  be  unscattered  signal  so  that  f(FOV)  =1).  With 
the  limited  data  available,  it  is  difficult  to  determine  an  analytical 
function  for  f  (FOV)  with  ar.y  great  precision.  However,  the  functions 
appear  to  be  roughly  logorithmic  with  FOV.  Assuming  this  to  be  the 
case,  our  empirical  formula  for  f(FOV)  (FOV  in  mr)  is  given  by 

f (FOV)  =  1  -  k  In (FOV)  (1  <  FOV  <  1500  mr) ,  (2.4-2) 


km  DATA  0.3  km  DATA 


NOISSIIAISNVdl 


Comparison  of  Transmission 
Data  for  0.3  km  and  1.6  km 
paths  (FOV  =  17°  Full  Angle) 


Experimental  Values  of  f (FOV) 
(Dashed  Lines  Bound  Upper  and 
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2.5  Relation  of  tJV  to  Visible  Extinction 

Of  crucial  importance  to  assessing  tne  expected  performance  of 
a  field  operated  communication  system  is  the  ability  to  make  relatively 
simple  optical  thickness  measurements.  In  the  visible  portion  of  the 
spectrum,  crude  estimates  of  optical  thickness  can  be  obtained  from 
observations  of  visual  range.  However,  to  use  such  observations  for 
a  uv  system,  it  is  necessary  to  relate  the  extinction  coefficients  for 
uv  and  visible  light.  To  this  end,  we  have  made  extensive  measurements 
of  middle  uv  (2537  ft)  versus  HeNe  laser  (6328  ft)  extinction  coeffi¬ 
cients  in  fog.  For  these  experiments,  the  two  sources  were  located  at 
a  distance  of  0.43  km  from  the  receiver  and  0.3  meters  apart.  A 
single  narrow  FOV  receiver  was  sequentially  aimed  at  one  or  the  other 
of  the  two  sources,  each  of  which  was  viewed  with  a  narrow  band  filter 
centered  at  the  appropriate  wavelength.  (Once  clear  weather  acquisi¬ 
tion  of  one  or  the  other  of  the  sources  was  obtained,  the  other  source 
could  be  acquired  easily,  even  in  heavy  fog.  Relatively  little  re¬ 
aiming  was  required  and  could  usually  be  accomplished  by  peaking  the 
signal  from  one  or  the  other  sources.) 


Fig.  2-15  is  a  plot  of  measured  uv  versus  HeNe  laser  extinction 

coefficients.  Observe  that  for  greater  than  8  km  1,  the  two 

coefficients  are  roughly  linearly  related.  Of  course,  since  the  two 

must  both  be  very  close  to  zero  in  clear  weather  (ignoring  ozone 

absorption  in  the  uv) ,  this  behavior  cannot  continue  for  smaller 

values  of  a  .  We  do  not  have  any  data  in  the  region 

HeNe 

aHer'e  <  8.  region  corresponds  to  the  onset  of  fog  which  usually 

occurs  rapidly  enough  that  it  is  very  difficult  to  make  reliable 
measurements).  In  the  region  8  <  a„  <  20,  a  least  squares  fit  to 
the  data  in  Fig  2-15  is 


a  -  -2.06  +  0.83  a  „  M  (2.5-1) 

uv  HeNe 

The  slope  of  0.83  relating  the  two  coefficients  is  somewhat 
disturbing  in  light  of  the  published  measurements  of  Baum  and  Dunkelman 
(2  8]  and  the  theoretical  predictions  of  Shettle  [29].  Shettle  pre¬ 
dicts  that  the  fog  extinction  coefficient  is  insensitive  to  wavelength 
in  the  region  from  the  middle  uv  through  the  visible.  Barur.  and 
Dunkelman  show  the  coefficient  to  be  slightly  decreasing  as 
wavelength  increases,  although  their  measurements  may  be  biased  toward 


to  Clear  Weather  Point) 
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higher  uv  coefficients  due  to  uncertainties  about  ozone  concentrations. 


In  both  cases,  the  slope  of  the 


a  vs . 
uv 


HeNe 


curve  is  ,  greater 


than  or  equal  to  one,  compared  to  a  slope  of  0.83  for  our  measurements. 

Note,  however,  that  the  error  bars  on  the  clear  weather  point 
(with  respect  to  which  all  the  data  are  plotted)  are  rather  large. 


spanning  a  distance  of  1.6  km 


-1 


for 


a 


HSNe 


and  +0.7  km  ^  for 


aHeNe  around  the  straight  line  fit.  Hence,  the  a 

-1  uv 

curve  may  be  shifted  to  the  left  by  2.3  km  or  down  by 
The  data  in  Fig.  2-15  are  thus  bounded  by  the  two  lines. 


auv  =  "°*15  +  °*83  aHeNe 


a 


uv 


-3.76  +  0.83  a 

HeNe 


(2.5-2) 


This  region  is  shown  in  Fig.  2-16  along  with  the  theoretical 
curve  (aHeNe  ~  auv)-  A  straight  line  with  a  slope  of  one  is  certainly 
possible  within  the  range  of  the  data  indicated.  However,  there  is 
still  a  small  absolute  shift  toward  smaller  values  of  a  uv. 

The  difference  between  our  measured  relationship  and  Shettle's 
predicted  one  may  reflect  the  difference  between  the  actual  parcicle 
size  distribution  we  encountered  and  his  specific  assumptions  about  the 
particle  size  distribution.  "lis  curves  show  that  there  is  some  wave¬ 
length  beyond  which  the  fog  extinction  coefficients  increase  slightly 
with  increasing  wavelength.  (For  his  choice  of  particle  sizes,  this 
critical  wavelength  is  above  1  ym. )  This  critical  wavelength  is 
most  probably  the  one  at  which  the  wave length  becomes  comparable  to  the 
mode  radius  of  the  particle  size  distribution.  If  the  distribution  were 
shifted  somewhat  toward  shorter  values  of  radii,  the  extinction 
coefficients  might  exhibit  the  observed  increase  of  1-2  km  ^  at  the  high 
end  of  the  visible  spectrum. 


(km 
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3.0  Comparison  with  Existing  Propagation  Theories 

In  this  section,  a  discussion  of  existing  channel  propagation 

theories  is  presented.  These  theories  include  the  extremes  of  single  scatter 

and  diffusion,  as  well  as  multiple  forward  scatter  theory1  [6,  16,  17]. 

Comparisons  between  the  predictions  of  these  theories  and  the  experimental 

data  presented  in  Section  2.0  is  the  primary  goal  of  this  section.  The 

comparisons  give  a  clear  idea  of  where  the  data  lies  in  the  spectrum  between 

diffusion  and  single#scatter ,  and  of  what  still  needs  to  be  understood. 

3.1  Results  from  Propagation  Theories 

In  this  section,  the  relevant  theoretical  expressions  are  presented, 
with  a  brief  discussion  of  the  appropriate  context  of  each  of  the  theories. 

3.1.1.  Single  Scatter  Theory 

Single  scatter  theory  adequately  describes  propagation  when  the 
probability  is  very  small  that  a  photon  will  be  scattered  more  than  once 
between  transmitter  and  receiver.  An  alternative  viewpoint  for  larger  optical 
thicknesses  is  that  single  scattering  theory  describes  the  first  scattering 
order  of  the  multiply  scattered  field. 

The  point  of  departure  for  the  present  discussion  is  the  theory  presented 
in  reference  [13 j,  based  on  the  work  in  reference  fi].  The  authors  present 
expressions  for  received  Intensity  as  a  function  of  time  for  an  impulsive 
source.  The  derivation  is  carried  out  in  prolate  spheroidal  coordinates,  with 
the  transmitter  and  receiver  located  at  the  foci  of  an  ellipsoid.  See  Fig. 
3-1.  This  approach  is  well  suited  to  the  problem  because  of  the  property  that 
the  sum  of  the  distances  between  the  foci  and  any  point  on  the  ellipsoid  is  a 
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constant.  Hence,  ellipsoids  are  equitemporal  surfaces;  i.e.,  computing  the 
total  single  scattered  radiation  reaching  the  receiver  at  any  one  time  only 

involves  an  integration  over  this  surface. 

The  total  received  energy  density  integrated  over  the  entire  scattering 

volume  is  given  by 


E 


r 


Vs 

.  v  4lTfit 


exp  [-ct(rt  +  rr)  ] 


p(0s)dV  , 


(3-D 


where 

I  Et-transmitter  energy  in  [joules]; 

1  Er=received  scattered  energy  density,  in  [joules/m2]; 

I 

t  o£  =extinction  coefficient; 

i  i. 

[  oto=scattering  coefficient; 

|  I  1 

I  r,=distance  from  transmitter  to  scattering  event; 

k  *  1 

[  rr=distance  from  scattering  event  to  receiver; 

i 


VC  w**"«b-* ' 
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JX.t=transmitter  solid  angle; 

A  =scattering  angle  at  the  scattering  event; 
p(  9  )=single  scatter  phase  function; 

3 

V=scattering  volume  as  determined  by  the  transmitter  beam 
pattern, 

The  phase  function  is  a  composite  of  two  functions: 


pRaylelghtf,s>  =  4  (1  +  008  V 


(3-2) 


wv  - 


i_Z_3L 


(1  +  g  -  2g  cos  8g) 


1.5 


(3-3) 


Thus, 


V.  pMie(9s>  +  «, 

p total  8s  =  - _ 


RayleiqhPRayle^h(9s> 

-  - - - ;  (3_4) 


a  +  a 
s  s 

Mie  Rayleigh 


where  <<  =  at  +  at  .  The  Mie  phase  function  is  the  commonly  used 

Mie  sRayleigh 

Henyey-Greenstein  phase  function  [18  ],  Here,  g  is  tho  average  cosine  of  the 
aerosol  scattering  function: 


9  ~  271  0  COS0sPMie(0s)sir‘e-ds* 


(3-5) 


The  single  scatter  theory  as  presented  in  reference  [18]  assumes  a  pulsed 
source.  For  comparison  of  the  single  scatter  results  with  cw  source 
measurements ,  the  received  pulse  is  simply  integrated  over  all  time. 
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3.1.2.  Diffusion  Theory 

Diffusion  theory  is  applicable  in  a  very  highly  multiple  scattering 
environment,  which  is  the  opposite  extreme  to  that  of  single  scatter.  The 
precise  point  at  which  atmospheric  conditions  are  such  that  diffusion  theory 
is  applicable  to  uv  propagation  is  still  an  open  question.  However,  an 
estimate  can  be  made  as  follows:  Kennedy  r.19 ]j  has  shown  that  for  isotropic 
scatter  with  large  single  scatter  albedo  (>0.9),  diffusion  theory  gives 
essentially  the  same  results  as  the  exact  transport  equation  solution,  for 
optical  thicknesses  between  1  and  100,  Bucher  [201  has  shown  that  many 
isotropic  scatter  results  can  be  applied  to  non-isotropic  scatter  situations 
as  long  as  the  diffusion  thickness,  is  equivalent  in  the  two  cases  and  is 
greater  than  3.  Here 


Td  =  t(1  -  g)  ,  (3-6) 

where  1  is  the  optical  thickness  of  the  medium  and  g  is  the  single  scatter 
phase  function's  average  cosine.  (See  Eq.  (3-5).)  For  an  average  cosine  of 
g=0.85  (typical  of  middle  uv  phase  functions  in  fog),  this  criterion  implies 
an  optical  thickness  greater  than  20. 

The  diffusion  theory  expressions  relevant  to  the  present  discussion  can 
be  derived  by  introducing  Pick's  diffusion  law  into  the  transport  equation 
[19,21  With  a  cw  isotropic  point  source  embedded  in  an  infinite  medium,  the 
received  irradiance  is  given  by  [19 j 
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^(L,8)  =  p  cos  0  exP  (-KL) 
A  (4tt)2DL 


n  j.  3D  ( 1  +  KL ) 

U  +  - - -  cos  Qj 


(3-7) 


where 

P^stransmitter  power; 
c=speed  of  light; 

L=source  to  receiver  distance  along  the  optical  axis; 

B  = angle  of  arrival  of  the  incoming  light  relative  to  the 
optical  axis; 

D=1/C3(rf-g«>Cg)]; 

<  =extinction  coefficient; 

«<  ^scattering  coefficient; 

K=[«-o<.  )/D]1/2. 

5 


Eq.  (3-7)  is  the  diffusion  expression  for  a  cw  source.  For  a  pulsed 
source,  the  time  dependent  solution  to  the  diffusion  regime  transport  equation 
yields  the  following  expressions  for  the  full  width  to  half  maximum  (FWHM) 
spread  of  the  received  pulse  [19]: 
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fl.ll} 


4cD 


FWHM  = 


(0.37  -  1.7) 


i.67 


if 

4  CD 
1/4 


c(cx  ) 

a 


37T 


aa  =  0 


a  L 
a 

4D 


^  1.5 


(3-8) 


Vi 

4D 


>  1.5 


where  4  =  <-  oi  is  the  absorption  coefficient  and  all  other  quantities  are  as 

H  3 

defined  following  Eq.  (3-7). 

In  subsection  3.2.3,  the  diffusion  theory  FWHM  pulse  spread  is  compared 
with  the  single  scatter  and  multiple  forward  scatter  results. 


3.1.3  Multiple  Forward  scatter  Theory 

The  narrow  angle  forward  scattering  approximation  has  been  receiving 
increasing  attention  in  atmospheric  scattering  contexts  [6,16,17,22  j  Although 
different  authors  introduce  it  in  a  variety  of  ways,  the  approximation 
essentially  consists  of  replacing  the  actual  single  scatter  phase  function 
with  one  that  has  only  a  very  narrow  forward  scatter  component.  All  the 
backscattered  and  wide-angle  scattered  light  is  assumed  to  be  completely  lost 
to  the  receiver,  and  hence  is  lumped  into  an  absorption  term. 

The  motivation  for  using  such  an  approximation  stems  from  the  persistent 

narrowness  of  measured  angular  spectra  and  lack  of  significant  multipath 

/< 

spreading  for  optical  thicknesses  less  than  20  il7,23,24,.  Bucher  has  observed 
[25]  that  the  dominance  of  the  multiple  forward  scattered  component  in  the 
received  field  may  exr’ain  the  better  than  expected  performance  of  many  simple 
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laser  communication  systems  in  bad  weather. 

The  specific  formulation  of  the  multiple  forward  scatter  theory  that  will 
be  used  here  is  the  Extended  Huygens-Fresnel  Principle  developed  by  Shapiro 
f 6  j  in  analogy  with  results  from  turbulence  theory  [26].  The  geometry  is 
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Fig.  3-2  Geometry  for  Extended  Huygens-  Fresnel  Principle 


shown  in  Fig.  3-2. 

The  discussion  in  reference  [  6  ]  is  exclusively  in  terms  of  a  Jaser 
source.  Here  the  basic  theoretical  result  is  applied  to  a  spatially 
incoherent  source  such  as  the  mercury  vapor  lamp  used  in  the  uv  experiments. 
For  the  conditions  encountered  in  the  experiments,  the  resulting  expression 
for  <P£j>.  the  average  power  detected  at  the  focal  plane,  is 


<Pd>  = 


<P_>  (FOV) 4 

K 


+  (FOV) 


2] 


(3-9) 


where  <P  >  is  the  average  power  arriving  at  the  receiver  entrance  pupil, 
R 
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where  <P  > 
R 


is  the  average  power  arriving  at  the  receiver  entrance  pupil* 


pt  2  "<xaL 
'V  ■  ~l  *  ps  e  ■ 

L 


FOV  is  the  full  angle  FOV  of  the  receiver,  and 


-1/2 


p0  -  (%L0ok  /3) 


(3-10) 


(3-11) 


is  the  channel  coherence  length.  6Q  is  the  width  of  the  narrow  forward  peak 
in  the  phase  function  (See  Fig.  3-3.),  and  k  =  2tt/A. 

Multipath  time  spread  in  the  multiple  forward  scatter  regime  is  given  by  [27] 


rms  width 


L  f  0-30 
C  0 


[1  +  2.250) 


1.5 


-  1] 


-  1 


(3-12  0 


Here  c  is  the  3peed  of  light  and 


0  =  (3J3  ' 

i3  the  rms  angular  beamwidth  of  the  received  signal.  is  the  phase 
function  width  introduced  above  and  *£  is  the  optical  thickness.  is  the 


single  scatter  albedo: 
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a. 


W0  = 


a  +  a 
a  s 


(3-14) 


One  final  point  should  be  made  about  utilizing  the  multiple  forward 
scatter  theory  with  real  phase  functions,  which  typically  have  significant 
wide-angle  scatter  and  backscatter  components.  Since,  as  seated  above,  it  is 
a  basic  premise  of  the  theory  that  these  non-forward-scatter  components  are 
never  collected  by  the  receiver,  they  nre  lumped  into  absorption.  Thus  the 
actual  scattering  and  absorption  coefficients  in  the  medium  are  scaled  as 
follows ; 


a 


I 

s 


aa  *  aa  +  ^  “  *)as 

where  <f>  is  the  forward  scatter  efficiency  tefined  by 


(3-15) 


4>  =  2tt 

ftp  is  the  effective  width  of  the  forward  scatter  peak.  (There  are  various 
ways  of  determining  The  method  used  here  is  that  used  in 

reference  [17)  The  quantities  ^  and  are  used  for  computations  with 
the  theory  instead  of  *  and  *  . 


P l 9 )  sin  ede 


(3-16) 


3.2  Application  to  UV  Experiments 

In  this  section,  the  theoretical  results  presented  above  will  be  used 
to  compute  the  expected  angular  spreading  for  the  conditions 
encountered  in  the  uv  field  experiments.  Comparisons  will  be  made  with  the 
measured  data.  The  three  theories  will  also  be  used  to  predict  the  multipath 
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time  spread  expected  under  these  conditions. 


As  a  preliminary  to  these  comparisons,  subsection  3.2.1  will  discuss  the 
applicable  phase  functions  needed  for  the  computations. 

3.2.1  Middle  UV  Phase  Functions  in  Fog 

The  accuracy  of  single  scatter  and  multiple  forward  scatter  computations 
depends  upon  the  use  of  an  appropriate  phase  function.  However,  in  most  case 
of  interest  it  is  parameters  such  as  the  width  and  zero-angle  peak  value  that 
are  most  important,  while  the  specific  shape  is  less  important.  Hence,  a 
variety  of  different  shapes  are  used  (e.g.,  Henyey-Greenstein  and  Gaussian). 
As  long  as  the  free  parameters  in  the  proposed  phase  function  are  adjusted  to 
embody  the  essential  features  of  the  physical  situation,  the  results  of  the 
theory  should  be  accurate. 

As  stated  in  Section  2.1,  no  particle  size  measurements  were  made  during 
the  uv  experiments.  Hence  it  is  necessary  to  rely  on  published  data  for 
comparisons  of  theory  with  our  experiments.  Because  of  the  relatively 
consistent  weather  conditions  encountered  during  the  experiments,  a  phase 
function  for  a  stable  maritime  fog  should  be  the  appropriate  one  to  use.  How¬ 
ever,  no  measurements  of  this  type  of  phase  function  in  the  middle  uv  have 
been  published.  Reference  (10  ]  contains  a  measured  fog  phase  function 

at  a  wavelength  of  250  nm,  but  it  is  for  a  radiational  (inland)  fog.  For  the 
purposes  of  initial  comparison,  the  phase  function  that  will  be  used  is  that 
presented  in  reference  [17J,  which  is  typical  for  0.53/Mm  propagation  in 
maritime  fogs  at  Point  Loaa,  San  Diego,  California.  (See  Fig.  3”3»)  Section 

3.2.2  will  include  a  discussion  of  how  the  theoretical  results  depend  on  the 
specific  phase  function  parameters  chosen. 

For  the  Point  Loma  phase  function,  reference  [17]  gives  the  following 
values  for  the  parameters  6q  and  S: 
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3.2.2.  Angular  spreading 

This  subsection  uses  the  theories  presented  in  Section  3.1  to  predict 
angular  spreading  behavior  for  the  experimental  conditions  encountered  in 


ANGLE  (DEGREES) 

Lubec  and  Eastport. 

Me; sured  and  theoretical  signal  versus  FOV  characteristics  are  plotted  in 
Fig.  3— **  For  two  values  of  X  •  The  large  FOV  data  is  in  two  groups:  the  lower 
group  of  data  represents  X=5  and  the  upper  group  X=9.  At  small  fields  of 
view,  the  data  was  not  very  sensitive  to  changes  in  X ;  the  entire  range  of 
values  observed  for  5<t<9  was  included  in  the  error  bars.  ,'This  is  consistent 
with  Figs.  2-8  through  2-10,  in  which  the  slopes  of  the  5.  10  and  15  ror  curves 
do  not  signif leant. !y  differ.)  Note  that  since  the  actual  signal  versus  FOV  mea¬ 
surements  included  the  unscattered  component,  this  has  been  subtracted  out  of  the 
data  in  the  plot.  Thus  the  data  points  represent  purely  scattered  signal.. 

It  is  immediately  evident  upon  inspecting  Fig.  3-M  that  both  the 
diffusion  component  and  t re  single  scatter  component  are  an  order  of  magnitude 
below  the  measured  data  for  fields  of  view  less  than  50  mr.  The  multiple 
forward  scatter  component  for  %=S  is  quite  close  to  the  measured  data  for 
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fields  of  view  less  than  15  mr,  and  both  the  £=5  and  curves  are  close  to 
the  data  at  larger  fields  of  view.  One  interestng  feature  of  these  curves  is 
that  while  the  multiple  forward  scatter  component  flattens  out  completely  for 
angles  greater  than  about  200  mr,  the  data  appears  to  be  increasing  in  that 
region.  Note,  however,  that  in  both  cases  the  diffusion  component  becomes 
significant  at  these  large  fields  of  view.  Thus  the  measured  large  FOV 
behavior  might  be  explained  by  a  change  in  dominance  from  the  multiple  forward 
scatter  component  to  the  diffusion  component. 

The  conclusion  that  diffusion  and  single  scatter  contribute  a  negligible 
amount  to  the  small  FOV  received  signal  is  quite  insensitive  to  specific 
assumptions  about  the  phase  function.  Fig,  3-5  show3  the  variation  in  the  two 
characteristics  as  the  parameter  g  is  varied.  The  optical  thickness  is  5. 

It  is  clear  from  these  curves  that  diffusion  depends  only  minimally  on  the 
value  of  g.  While  the  single  scatter  variations  are  more  pronounced,  the  most 
extreme  case  is  still  below  the  unscattered  component  for  fields  of  view  less 
than  25  mr. 

The  variation  in  the  multiple  forward  scatter  component  is  somewhat  more 
complicated  because  of  the  dependence  on  the  two  parameters  and  ^ . 
Re-writing  Eq.  (3-9)  in  terms  of  these  parameters  and  dividing  by  the 
unscattered  component,  shows  this  dependence  explicitly: 


<Pd> 


2  T^u0 
(FOV)  e  u 


P  2  2 

unscattered  5.4890<J>t  +  (FOV) 


(3-17) 


Here  FOV  is  the  full  angle  field  of  view  and  is  the  scattering  albedo. 
Observe  that  the  exclusive  role  of  is  to  determine  the  breakpoint  in  the 
curve,  while  f>  affects  both  the  breakpoint  (to  a  lesser  extent  than  dQ)  and 
the  final  value  ( exponentially) . 

Figures  3-f  and  3-7  show  multiple  forward  scatter  signal  versus  FOV 
characteristic  for  various  values  of  and  <f>  .  Note  from  Fig.  3-6  that  0 


Fig.  3-5:  Influence  of  Average  Cosine 
on  Single  Scatter  and 
Diffusion  Results 
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does  not  significantly  affect  the  small  FOV  behavior  of  the  curve,  although  it 
changes  the  final  value  by  almost  an  order  of  magnitude,  From  Fig.  3-7, 
however,  it  is  clear  that  drastically  affects  the  small  angle  signal  by 
changing  the  breakpoint  of  the  curve.  Because  of  this  critical  dependence  on 

it  is  difficult  to  say  that  the  multiple  forward  scatter  theory  matches 
our  experimental  data.  By  appropriate  choice  of  and  <f,  the  theory  can  be 
made  to  match  a  wide  range  of  data.  Thus,  an  accurate,  independent 
determination  of  these  parameters  would  be  required  to  state  with  certainty  that 
the  theory  is  a  good  fit  to  the  data. 

On  the  other  hand,  the  narrow  forward  angle  peak  of  the  fog  phase 
function  coupled  with  the  observed  narrowness  of  the  angular  spectrum, 
strongly  recommends  the  multiple  forward  scatter  theory  or  something  similar. 
And,  of  course,  of  the  th-ee  theories  discussed,  the  multiple  forward  scatter 
theory  is  the  only  one  that  predicts  anything  like  the  observed  behavior.  It 
is,  in  fact,  remarkable  that  with  all  the  uncertainties  involved  in  using  the 
Point  Loma  phase  function,  the  signal  versus  FOV  characteristic  for  £Q=11.3 
mr,  0.57  and  t=9  is  reasonably  close  to  the  data  for  most  fields  of  view 
observed . 

One  final  remark  is  appropriate  nere  before  proceeding  with  the 
discussion  of  multipath  time  spread.  If  the  multiple  forward  scatter 
component  proves  to  be  the  dominant  one  for  narrow  receiver  fields  of  view, 
this  has  certain  implications  for  the  difference  between  collimated  and  wide 
angle  sources.  In  particular,  it  implies  that  the  character  of  the  received 
field  is  not  significantly  different  for  the  two  types  of  sources,  since  only 
light  which  remains  confined  to  a  small  cone  around  the  optical  axis  has  a 
chance  of  reaching  the  receiver. 
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3.2.3  Multipath  Time  Spread  Estimates 

This  section  provides  a  brief  discussion  of  multipath  time  spread  and 
estimates  of  time  spread  are  made  for  the  experimental  conditions  encountered 
in  the  uv  field  trips. 

In  general,  measured  multipath  time  spread  depends  on  receiver  FOV  as  well 
as  channel  conditions.  This  is  a  consequence  of  the  larger  percentage  of  longer 
photon  paths  among  those  which  make  up  the  off-axis  signal.  This  angular 
dependence  varies  among  the  different  components  of  the  field  discussed  above. 

On  one  extreme  is  the  very  explicit  relationship  between  FOV  and  time  spread  in 
the  single  scatter  component.  For  the  coaxial  geometry  illustrated  in  Fig  3-1, 
the  intersection  of  the  transmitter  and  receiver  cone  angles  sets  an  absolute 
upper  limit  to  the  time  spread.  On  the  other  extreme  is  the  diffusion  component, 
which  has  no  angular  dependence. 

Between  these  extremes  is  the  multiple  forward  scatter  component.  Observe, 
however,  that  the  time  spread  expression  derived  by  Stotts  [ Eq . ( 3—12) ]  does  not 
include  any  FOV  dependence.  Upon  inspection  of  his  derivation  [27],  it  is 
evident  that  Eq.  (3-l2)applies  to  the  situation  in  which  the  receiver  FOV  is 
large  enough  to  collect  all  of  the  multiple  forward  scatter  component.  Hence, 
his  expression  gives  an  upper  bound  to  the  multipath  time  spread  attributable 
to  this  component. 

Table  3-2  shows  computed  values  of  time  spread  using  the  three  theories., 

Bucher's  simulation  formula  [20]  is  also  included.  The  physical  conditions 
assumed  for  these  estimates  are: 

L=.31  km 
g=0. 85 
(J0=0.96 

4-0.51 

£>0=11.3  mr 

FCV  (full  angle)= 100  mr 


The  FOV  is  taken  to  be  large  enough  so  that  the  majority  of  the  multiple 


forward  scattered  light  is  collected, 

Table  3-2.  TIME  SPREAD  ESTIMATES  FOR  0.31  KM  PATH  (SEC) 


X 

Single 

Scatter 

(FVJHM) 

MFS 

(ras) 

Bucher 

(rms) 

Diffusion 

(FWIifl) 

5 

5x10"n 

1 .31x10**® 

5. 14x1C~7 

2.5-11. 9x1 0~7 

7 

5X10"11 

1.30x10"8 

6.77x10“7 

3. 6-16. 5x1 0-7 

9 

5x10-11 

1.31x10"8 

8.3x10"7 

4.7-21 . 4x1 0~7 

Observe  from  the  table  that  even  in  the  worst  case  of  diffusion,  the 
multipath  spread  does  not  significantly  exceed  2/isec.  This  is  consistent 
with  the  lack  of  time  spreading  observed  with  the  i-2 ytosec  pulsed  source  in 
the  Eastport  experiments  for  optical  thicknesses  less  than  10. 
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4.0  Transport  equation  simplifications 

One  major  purpose  of  the  uv  experiments  has  been  to  obtain 
enough  information  about  the  propagation  parameters  so  that  the 
transport  equation  can  be  simplified.  In  this  section,  we  discuss 
one  particular  simplification,  which  results  from  the  observed 
insensitivity  of  the  scattered  signal  component  to  spatial  varia¬ 
tions  in  the  medium's  extinction  coefficient  (see  Section  2.4). 

(Another  simplification,  based  on  the  relative  narrowness  of  the 
measured  angular  spectrum,  is  discussed  elsewhere  [23 j.)  Specifi¬ 

cally,  it  will  be  shown  that  under  certain  conditions  only  the  homo¬ 
geneous  transport  equation  needs  to  be  solved.  The  basic  argument 
and  result  is  presented  here,  but  no  detailed  solutions  are 
pursued. 

Consider  the  generally  spatially  inhomogeneous  linear  transport 
equation  [7,8] 

(a  (r)  +  ^  “  +  ft  •  v)  f  (r,  ft,  t) 

«  as  (r)  fdH  p(ft*ft  )  f(r,  ft  ,  t) , 

where  £(r,  ft ,  t)  is  the  probability  density  of  time  t  for  a  photon 

to  be  at  point  r,  going  in  direction  ft; 

a  (r )  is  the  extinction  coefficient  at  point  r  'in  the. 

_  medium; 

ag (r)  is  the  scattering  coefficient  at  poin_  r  in  the 

medium; 

_  —  * 

p(ft*ft  )  is  the  normalized  single  scattering  phase  function; 

c  is  the  speed  of  light; 

7  is  the  gradient  operator. 

Observe  that  the  angular  characteristics  of  the  phase  function  are 
assumed  to  be  the  same  throughout  the  medium;  the  inhomogeneity  is 
lumped  into  che  scattering  coefficient  c»s  (r)  .  Both  a(’r)  and  as(r) 
are  random  processes. 

To  simplify  the  subsequent  discussion,  it  will  be  helpful  to 
consider  the  integral  equation  corresponding  to  Eq(4-1)  in  uhe 
Laplace  domain.  Taking  the  one-sided  Laplace  transform  of  Eq  (4-1) 
and  integrating  with  respect  to  the  pathlength  variable  in  the  direction 
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of  photon  travel  gives 


F(r,  ft,  s) 


/ d  ag(r  -  vft)  / d  'p(n-fl')F(r  -  vft,  ft',  s) 

o  3 

(4-2) 


•exp  [-  J  y  (r  -  v  '  ft)  dv  '  ] 
0 


17  ___  v_ 

+  -  j  dy  Fq  (r  -vft,  $7)  exp  [  —  /  y  (r  -  v'ft)dv'], 
0  o 

where 


CO 

F(r,  ft,  s)  =  /  ds  f(r,  ft,  t)  e  st 

(4-3) 

0 

y (r)  =  a (r)  +  s/c, 

(4-4) 

and 


F0(r,  ft)  =  f  (r,  ft,  t  »  0)  ,  (4-5) 

The  second  term  in  Eq  (4-2)  is  a  source  term  introduced  by  the 
initial  value  property  of  the  one-sided  Laplace  transform.  It 
corresponds  to  the  unscattered  component  of  the  received  field, 
and  is  denoted  by  Fu  (r,  R,  s) 

The  insensitivity  property  of  the  scattered  field  is  now 
introduced  in  the  following  manner:  Eq.  (4-2)  is  converted  to 
an  integral  equation  for  the  scattered  part  of  F(r,  Q,  s)  , 
denoted  by  Fg  (r  ,17,  s).  This  is  accomplished  by  substituting 

F(r,  fl,  s)  =  F  ( r ,  ft,  s)  +  F  (r,  ft,  s)  (4-6) 

u  s 

on  both  sides  of  the  equation  and  cancelling  common  terms.  The 
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result  is: 

00 

F  (r,  IT,  s)  =/  dv  a  (r  -  vST)  /  dft'p($7  •  s2')F  (r  -  vii7,  ft1  ,  s) 

S  o  s  s 

V 

exp[-f  y (r  -  V  ft)  dv' ] 

0 

oo 

+  /  dva  (r  -  v  ft)  Jdft’  p(ft*ft')F  (r  -  vft,  ft'  ,  s) 
s  u 

0 

V  _ 

exp [-/  a(r  -  v'ft)dv'].  (4-7) 

c 

The  first  term  in  Eq  (4-7)  is  the  multiple  scatter  portion  of 
F  (r,  ft,  s)  and  the  second  term  is  the  single  scatter  portion. 

It  can  be  shown  that  under  certain  conditions,  the  second  term 
in  Eq  (4-7)  has  the  same  dependence  on  the  variations  in  a(r)  as 
does  the  unscattered  component  itself.  This  implies  that  this  term 
must  be  small  compared  to  the  first  term  since  the  left  hand  side  of 
Eq  (4-7)  must  be  insensitive  to  these  variations.  Furthermore,  the 
first  term  must  also  be  insensitive  to  these  variations,  so  that  it 
must  behave  nearly  as  if  its  randomly  space-varying  extinction  and 
scattering  coefficients  were  replaced  by  their  means.  Hence,  Eq.  (4-7) 
can  be  written  as 

vft,  ft',  s '  e  ~J<a>v 
(4-8) 

(Two  basic  assumptions  are  implicit  in  this  argument:  first,  that  there 
is  no  exotic  form  of  cancellation  occurring  between  the  single  and 
multiple  scattering  terms;  and  second,  that  the  first  term  is  not  the 
small  term,  i.e.,  this  is  not  simply  a  single  scatter  problem.) 

Sufficient  conditions  for  Eq  (4-8)  to  be  valid  are: 

^  is  not  vaery  far  off  the  optical  axis  (Q  •  —  =  1)  ; 

|r  j 

2)  p  ( £2 •  52  ’ )  has  a  pronounced  forward  peak  (p  (90°  )  <<p  (0°)  )  . 


Fs(r' 


s)  = 


<a  > 
s 


/  dv 
o 


/d£2’p 


fl*  )F 
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The  second  condition  is  not  very  restrictive,  and  is  satisfied  by  all 
for  phase  functions  (see  Fig.  3-3} .  The  first  condition  is  quite 
restrictive,  but  we  believe  it  can  be  relaxed.  On  general  physical 
grounds,  it  can  be  argued  that  the  further  off  axis  the  receiver 
looks ,  the  more  the  light  it  collects  will  be  dominated  by  multiple 

scatter. 

Furthermore,  the  experimental  evidence  (see  Fig.  3-4  and  3-5) 
is  consistent  with  this  argument.  Hence  it  is  reasonable  to  assume 
that  Eq  (4-8)  is  valid  for  any  ft. 

If  the  spatial  Fourier  transform 

^  _  _  _  _  -jrr 

F  (u),  ft,  s)  *  [dr  F  (r,  ft,  s)e  (4-9) 

s  J  s 

V 

is  introduced  into  Eq  (4-8) ,  the  integral  over  v  can  be  carried  out, 
resulting  in 

«x  >  _  ^ 

F  (w,  ft,  s)  =  - JdftVfft-fi’ )F  (w,  ,  s).  (4-10) 

s  <Y>  +  ju)  •  ft  s 

Eq.  (4-10)  is  the  desired  result.  It  is  a  singular  eigenfunction 

equation  for  1*  (w,  ft,  s) .  Although  time  has  not  been  available 

s 

to  explore  solutions  of  Eq(4-10),  it  is  essentially  identical  to  one 
considered  by  Case  and  Zweifel  [7] .  The  lack  of  the  source  term  is 
a  significant  simplification  over  the  original  transport  equation,  and 
may  afford  closed  form  solutions  for  physical  conditions  under  which 
the  assumptions  leading  to  Eq(4-10)  apply. 
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